Motion of grains, droplets, and bubbles in fluid-filled nano-pores 
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Molecular dynamics studies of nono-sized rigid grains, droplets and bubbles in nano-sized pores 
indicate that the drag force may have a hydrodynamic form if the moving object is dense and small 
compared to the pore diameter. Otherwise, the behavior is non- hydrodynamic. The terminal speed 
is insensitive to whether the falling droplet is made of liquid or a solid. The velocity profiles within 
droplets and bubbles that move in the pore are usually non-parabolic and distinct from those corre- 
sponding to individual fluids. The density profiles indicate motional shape distortion of the moving 
objects. 



I. INTRODUCTION. 

Classical problems of motion of objects through flu- 
ids continue to be of interest because of their multiple 
applications in natural sciences and because of their fre- 
quent relation to fundamental issues. There is a vari- 
ety of possible physical situations. The moving objects 
may may be rigid or fluid. They may also span many 
length scales: from macroscopic, such as snowflakes 
or sediments B, to micron-sized, such as water droplets 
and aerosols in clouds |||| . A new frontier arises in the 
context of micro- and nano-scale machinery |^|| which 
sometimes involves flows through narrow pores. Two- 
phase flows in micro- and nano-pores would require un- 
derstanding of flows and interactions between nano-sized 
droplets. One recent example of this situation is provided 
by the surface force apparatus studies of capillary con- 
densation in a nanoscale pore J7|. Droplets, bubbles, and 
grains of various sizes may move at various characteristic 
speeds. At large Reynolds numbers, smooth trajecto- 
ries may become replaced by various classes of complex 
trajectories, such as observed in motion of falling disks 
H or bubbles J|,|l(|. Another complexity may appear 
due to the droplet /bubble shape deformation, possibly 
combined with splitting and stacking PJPJI or due to a 
bubble collapse. 

In this paper, we focus on issues arising at the molec- 
ular level aspects of such problems by considering nano- 
sized objects moving in larger but still nano-sized pores 
at low Reynolds numbers. Our interest is in probing 
physics of nano-scale drag phenomena in simple molec- 
ular models and in assessing the validity of continuum 
physics concepts at this length scale. Our main purpose 
is to compare motion of droplets and bubbles to that 
of the solid spheres. The spheres are built of fluid-like 
atoms that are tethered to amorphously located centers. 
Our studies are based on the molecular dynamics (MD) 
approach 1 13- 171] which allows us to monitor motion of 



individual atoms but its use is restricted to nano-second 
time scales. Such atomic models should apply literar- 
ily only to the microscopic-scale systems but they may 
also be considered as toy models of larger droplets. Note 



also, that rupture or coalescence of droplets ||16| , |l8| -p0|, 
no matter how big, necessarily comes to a stage in which 
the interface morphology involves subcontinuum physics 
of nano-scale tendrils. 
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FIG. 1. A snapshot of solid sphere moving through a fluid 
filled cylindrical pore. Only the solid and wall atoms are 
shown for clarity. The solid consists of 454 tethered atoms 
that are tethered to amorphously placed centers. There are 
2212 pore wall atoms which are fixed rigidly. The space be- 
tween the solid and the pore is filled with 2281 fluid atoms 
(the case of fluid B). The left-hand figure represents a stere- 
ographic projection whereas the right-hand figure represents 
a projection from the top. 

A natural way to study motion of an object in a fluid is 
to apply a driving force and determine a drag force, 
that brings the motion to a stationary state. The simu- 
lations on hard sphere systems ]^l|j22|j2^] have indicated 
validity of Stokes's law [E3[ , 



Fd = 6nRr/v, 



(1) 



down to atom sized objects, when average quantities are 
monitored. In the Stokes law above, R, rj, and v denote, 
respectively, the radius of the moving sphere, the viscos- 
ity of the surrounding fluid, and the velocity of the body 



relative to an undisturbed fluid. However, the continuum 
mechanics predictions pE) | break down and lead to diver- 
gent forces when one considers an approach of a sphere 
to a fixed wall at a constant velocity. The MD results 
p4| . on the other hand, give evidence for a scenario in 
which the fluid atoms escape the region that is squeezed 
between the sphere and the wall, preventing any buildup 
in the force. 

This paper provides a pilot MD study of the drag force 
phenomena that relate to the motion of droplets and bub- 
bles. This study is aimed at making a comparison to a 
similar MD analysis of solid nanoscale objects performed 
by Vergeles et al. |24|]. We demonstrate that the drag 
force on droplets is similar to that on solid objects but the 
droplets undergo motional distortion. The MD results 
for the motion of droplets essentially agree with hydro- 
dynamic predictions but this is not so in the case of the 
bubbles. The bubbles that we study have densities which 
are so low that they correspond to the crossover region 
between the sub-continuum and continuum physics |26[] 
and a non-hydrodynamic behavior, at such small length 
scales, is not surprising. Another interesting finding is 
that the surface tension effects at nano-scale confinement 
generally yield velocity profiles such that the velocity in 
the single fluid regions is quite distinct from that corre- 
sponding to one fluid flows. 

We start our analysis in Section 2 where construction 
of the nanoscale pore is described and properties of the 
fluid that fills it are established. The overall geometry 
is illustrated in figure |l| and the z-direction corresponds 
to the pore axis. The primary role of the pore is to pin 
the fluid to its walls and thus to provide a stationary 
environment in which motion of falling objects can be 
studied. Our objective in this paper is focused on com- 
paring various drag force phenomena at nanoscale length 
scale, as discussed within a simple molecular model of a 
generic nature, and not the boundary phenomena aris- 
ing near the confining wall. We take then the wall to 
be built of Lennard-Joncs atoms fixed rigidly on an fee 
lattice and adjust parameters that minimimize layering 
in the density profile near the wall. The reason is that 
it seems proper to consider first the simplest kind of the 
drag phenomena in which special wall-related effects are 
missing. Studies of near-the wall physics and of excita- 
tions in the solid due to a flow would require considering 
a more realistic model, such as used in studies of fluid 
flows in carbon nanotubes by Tuzun et al. p7j]. 

We consider two kinds of fluids: one at a dense liquid 
density (A) and another at a dense gas density (B). Both 
are described by the Lennard-Jones potential 



ULj(r) = 4e 



(2) 



where e and a are the units of energy and distance respec- 
tively (for krypton (7=0.357 nm and e/fc B =201.9 K). The 
parameter d a ^ is f for a = /3=A and 0.5 for a = f3=B. 
This choice makes B more gas like. We adopt the cut- 
off of 2.5a in the potential. Our choice of parameters 



of the system, such as temperature, densities, and those 
related to the geometry and potentials of interactions is 
ultimately geared towards generating a situation, by trial 
and error, in which the densities of fluids A and B are 
roughly the same in both the droplet and bubble cases. 
Thus the droplet made of fluid A that moves in fluid B 
is approximately a symmetric image - density wise - of 
the bubble of fluid B that moves in fluid A. 

In Section 3, we discuss gravity driven motion down the 
z-axis of a spherically shaped atomic solid. We present 
results on the density and velocity fields of the surround- 
ing fluid and confirm the validity of Stokes's law. In 
Section 4, we discuss static properties of droplets of liq- 
uid A surrounded by fluid B and those of B-fluid bubbles 
immersed in liquid A. We consider the immiscible case in 
which d a,b=0. In Section 5, we give illustrative examples 
of diffusional merger of the droplets in an environment 
unconstrained by any pore walls. In Sections 6 and 7, 
we return to the pore geometry and consider motion of a 
single droplet and a single bubble under the influence of a 
gravitational force, g. We consider two cases: either the 
g is applied only to the molecules of the droplet /bubble 
(Section 6), or it is applied to all of the molecules of the 
system (Section 7). In a true nano-pore, the latter is 
more realistic. On the other hand, as a toy model, the 
former is easier to understand because this case mimics 
the situation in which the object moves in a stationary 
medium and in which only the immediate neighborhood 
of the object is disturbed. In other words, the case of the 
force acting only on the object, corresponds conceptually 
to a motion in a much larger sized fluid than available in 
MD simulations. 



II. THE SINGLE FLUID FILLED NANO-PORE 

The geometry of the pore is shown in Figure [j]. The 
pore is constructed by first generating an fcc-like lattice 
in which subsequent sites in the x, y, and z-directions dif- 
fer by a. The wall sites are obtained by cutting out an 
annulus of the inner and outer radii of 12a and 15cr re- 
spectively. The pore's axial length is 24a and the periodic 
boundary conditions are adopted along the z-direction. 
The primary purpose of introducing the walls is to pin the 
fluid at the walls and allow for dissipation. The periodic 
boundary conditions in the x and y directions would not 
lead to establishment of the terminal speed. For simplic- 
ity, we keep the pore wall atoms frozen at their crystalline 



positions. 




FIG. 2. The top panel shows the density profiles for fluids 
A and B, as indicated, obtained in a vertical cut through 
the pore. The cut is at y—0 and the data are averaged 
over the axial z direction. The bottom panel shows ve- 
locity profiles for fluids A and B under gravity driven flow 
with g = O.Ole/a. The corresponding central densities are 
p = 0.2/cr 3 and p = 0.75/cr 3 . 

As is known from the MD studies of channel flows 
pq [30f occuring between two parallel plates, the walls 
may induce a layered structure in the density profile on 
the scale of several a in the neighborhood of the walls. 
The magnitude of this effect depends on the strength of 
the interaction between the fluid and the wall and on the 
fluid density. We adopt the Lennard-Jones interaction for 
the wall-fluid interaction with e replaced by e w j = 0.7e. 
Our intention here is to provide some pinning (with a 
small slip length) and, at the same time, to minimize ex- 
istence of the layered structure so that the density profile 
is as flat as possible. This particular values of e w j was 
obtained by trial and error. Figure |^ shows the result- 
ing density profiles for fluids A and B, obtained with 
6554 and 2646 molecules respectively. The densities at 
the central region are correspondingly 0.75 and 0.20 a~ 3 . 
For fluid B, there is essentially one layer next to the wall 
and just a remnance of the second atomic layer beyond 
which the profile is flat. Larger values of e w j would 
result in establishment of several layers. For liquid A, 
the intermolecular attraction is stronger than that due 
to the wall atoms which results in a depletion layer near 
the wall but the flat region in the center is large enough 



to accomodate a relatively large sized bubble of fluid B, 
as will be discussed later. The slip length for fluid B is 
smaller than for liquid A even though B is much more 
rarefied |3(J. This is due to the different values of d a ,p 
in equation 2 but the same wall-fluid interactions for the 
two fluids. 

A characteristic time scale r = o^m/e) 1 / 2 , where m is 
the mass of a fluid molecule (the masses of molecules, in 
fluids A and B are idebtical) corresponds to the period of 
oscillations in the Lennard-Jones potential minimum of 
fluid A. The integration step was 0.005r and the starting 
velocity distribution is Maxwellian. Our thermostatting 
procedure was based on the Langevin noise [fjlf which 
balances the frictional dissipation. The equation of mo- 
tion for the x-coordinate of a molecule reads then 



mx 



F c -(i + r 



(3) 



where F c is the force due to other molecules, £ = 
0.005m/r is the coefficient of friction, and T is the Gaus- 
sian uncorrelated random force such that 



(T(t)T(t')) = 2(T5{t - t') 



(4) 



A similar equation holds for the other coordinates but 
no Langevin noise was applied in the z-direction in order 
not to affect the systematic flows due to gravity, once 
it is switched on. The equations of motion are solved 
by means of the fifth order predictor-corrector scheme 
||14f . Throughout the paper, we report on the calculations 
done at the temperature T of 0.71 e/ks at which there 
is a sufficient stability of the droplets of the size studied 
here. Another reason for this choice of T is that this is 
one of the several values chosen in the classic MD studies 
of droplets by Thompson et. al. p2fl . 

The viscosity of fluids A and B can be obtained from 
the velocity profiles, as shown if figure ^, by fitting to 
p3[ to the parabolic form 



P9 (U 2 2n 



(5) 



Here, Rq is the effective radius of the cylinder, p denotes 
fluid density, and r is the polar coordinate. Our values 
of n for fluids A and B are 2.22 ± 0.12 and 0.23 ± 0.04 

' (XT 

respectively. 



III. MOTION OF SPHERICAL SOLID GRAINS 

We begin our analysis by considering a gravity driven 
motion of a spherically shaped solid object in fluid B. Our 
model of the falling object is an amorphous variant of a 
model considered by Vergeles et. al. [Q. We first take 
454 points and place them randomly within a spherical 
volume with the density of 0.8/cr 3 . The corresponding 
radius is 5.1cr. These points are adopted as tethering 
centers for molecules of the same mass as that of the 
fluid molecules. The amorphicity is adopted for a more 
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direct comparison to a liquid droplet which is amorphous 
but lacks the stiffness of the solid The molecules of the 
solid are attracted to their tethering centers by a strong 
harmonic force with the spring constant of 400 e/er 2 . 

The solid atoms in the bulk of the sphere just stay at 
the tethering centers but those near the surface, in ad- 
dition to the elastic forces, experience interactions, u s j, 
with the molecules of the fluid. We chose u s j to be as 
between the fluid atoms, i.e. of the Lennard- Jones type: 



u sf (r) = 4e 



- d a,fi ~ 



(6) 



We have considered 3 values of d s j: 0, 0.5, and 1. The 
first choice, of a pure repulsion, mimics the interactions 
between the A and B fluid molecules to allow comparisons 
with the droplet and this is the choice we focus on more. 
The other choices allow the sphere to drag more fluid and 
perturb the environment stronger. 

We determine the motion of the solid atoms in two 
stages. First, we calculate what are the net force and 
torque with which the surrounding fluids acts on the 
whole solid. The tethering centers are then translated 
and rotated, around the center of mass, according to the 
values of these quantities. (We have found no systematic 
rotation in our studies). In the second stage, we reevalu- 
ate the forces on solid's atoms (the elastic contributions 
come from tethers whose anchores have been moved) and 
then evolve the positions of the atoms according to the 
standard MD scheme. Gravitational forces, if any, are 
applied after an equilibration period, of order 300r. 

Figure || shows the velocity of the center of mass as a 
function of time. The top panel is for d s j = and the 
bottom panel for two larger values of d s j . We observe the 
phenomenon of saturation - the center of mass velocity 
saturates at a terminal velocity Vt- The saturation level 
is a linear function of g. and it decreases when d s f is 
increased. The time scale to reach the saturation is of 
order 200r. From equation |l| for the grain we have found 
value of terminal velocity to be equal: v t = 0.203t/(t (g — 
O.Ole/cr) and v' t = 1.019r/a (g = 0.05e/cr). These values 
are in good agreement with ones obtained numerically 
when gravitational force is acting only on the atoms of 
the solid. For other case when his force being applied to 
all atoms the equation [I] is not valid. 

When the gravitational force is applied only to the 
atoms of the solid, the velocity field in the fluid is affected 
almost exclusively in the central region of the width that 
coincides with the sphere's diameter. If the force is ap- 
plied to all atoms, the terminal speeds are larger since 
the fluid itself participates in the motion - this becomes 
a Poiseuille-like flow. 
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FIG. 3. Evolution of the velocity of the center of mass of 
the spherical solid in fluid B for the two different values of 
g as indicated in the figure. The solid lines correspond to 
the gravitational force acting only on the atoms of the solid 
whereas the dotted lines correspond to this force being applied 
to all atoms. The upper panel is for d s ,/=0.5 and the terminal 
velocities are, top to bottom, 2.148, 1.297, 0.441, and 0.266 
a jr. The lower panel is for g = O.Ole/u and it compares two 
values of d s> f. 0.5 and 1. The terminal speeds are 0.254 and 
0.194cr/r respectively. 
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FIG. 4. The velocity profiles for the rigid spherical solid in 
a reference frame that moves with the grain and at the center 
of the grain. The solid lines correspond to the gravitational 
force acting only on the atoms of the solid whereas the dotted 
lines correspond to this force being applied to all atoms. The 
values of g, in units of e/o, are indicated. 



In Figure |j profiles of velocity of spherical solid are 
shown. These profiles have been obtained in the moving 



system of coordinates connected with the center of mass 
of spherical solid. The atoms in the central region move 
rigidly. 



IV. STATIC DROPLETS AND BUBBLES 

There have been a number of MD studies of droplets. 
Among the more recent entries there are simulations of 
droplets that form during a rapid expansion of a liq- 
uid [Q and the analysis of the terraced spreading of 
a 2000-molecule droplet that lands on an atomic surface 
p5| . A comprehensive study of liquid droplets in equilib- 
rium with their own vapor has been done by Thompson 
et. al. |52). They have considered between 41 to 2004 
Lennard- Jones molecules of a single kind. Still larger 
droplets have been simulated by Powles et. al. p6| , p7t - 
The studies by Thompson et. al. focus on the quantities 
related to the phenomenon of the surface tension. The 
authors have provided a molecular level validation for 
a) Laplace's equation for the pressure difference between 
the inside and outside of a droplet, b) Kelvin's equation 
for the radius-dependence of the vapor pressure, and c) 
Tolman's equation for the effect of curvature on the coef- 
ficient of the surface tension. The profiles of the density 
obtained on radial crossing the gas-liquid interface have 
been found to be well represented by the function: 

P(r) = \{pi + Pv) -\{Pi- Pv)tanh(2(r - R)/D s ), (7) 

where D s is a measure of the thickness of the interface, 
R is an estimate of the droplet radius, pi is the density 
of the liquid in the center of the droplet, and p v denotes 
the density of the vapor far away from the interface. D s 
depends on the gas-liquid interaction and it is usually of 
order several a. 

One theoretical device that was used by Thompson 
et. al. (3^] was to adopt the geometry of a spherical 
container that would scatter the molecules back to the 
system if they were about to leave it. When we studied 
droplets in a 'vacuum' with the periodic boundary con- 
ditions in all directions, we have found that at each value 
of T and of pi there is a threshold number of molecules 
below which the droplet evaporates and occupies all of 
the periodic volume. Above this threshold number the 
droplet remains stable and is surrounded by a few vapor 
molecules. 

We have checked that at T—0.7le/kB a 454-big droplet 
of atoms A, with an initial central density of about 0.8/a 3 
is stable on its own - without any spherical container. 
The stability properties are even more assured when such 
a droplet is immersed in the immiscible fluid B and this 
size of the droplet has been adopted for further studies. 
The confining pore walls also add to the stability. Figure 
H shows the equilibrated density profile when the droplet 
is placed in the pore filled with 2281 molecules of fluid 
B. We observe that the droplet is well defined and is not 



immediately interacting with the pore walls. Its radius 
is about the same as of the spherical solid studied in 
the previous section - about 5.1 er. The profile at the 
interface is consistent with the tanh form given by eq. [?]. 
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FIG. 5. The profile of the total density corresponding to a) 
a droplet of liquid A immersed in fluid B (the solid line), b) a 
bubble of fluid B immersed in liquid A (the dotted line), and 
c) the spherical solid, of the same number of particles as the 
droplet, immersed in fluid B. The molecules occupy the pore 
space and the axial average has been performed - the profile 
is plotted vs. the radius. The droplet data show an enhance- 
ment in the density next to the wall which corresponds to the 
first monolayer peak indicated in Figure [| without the axial 
averaging. 

In order to generate a stable bubble, we surround 300 
of the B-type molecules by 6145 A-type molecules. The 
equilibrated density profile is shown in Figure |. The 
corresponding radius is about 7.4c Thus the bubble is 
several ct's away from the depletion zone. The density 
(and pressure) within the bubble is essentially uniform 
and it is around 0.2tr 3 - similar to that of the fluid in 
which the droplet and the grain move. 



V. COALESCENCE OF DROPLETS 

Before we continue with our analysis of the drag force 
phenomena in the pore geometry we pause for some 
discussion of coalescence of two initially independent 
droplets. Continuum mechanics provides no mechanism 
for coalescence and yet, two droplets placed sufficiently 
close to each other may combine simply through atomic 
diffusion. 

In order to show this, we first consider a volume of size 
25 x 25 x 25 with the periodic boundary conditions. We 
form a single A-type droplet of 454 molecules in the cen- 
ter of the volume and surround it by 3751 molecules of 
type B (with the density of 0.2 cr -3 ). After an equilibra- 
tion period of 300r, we double the volume of the system 
by placing its mirror replica in the z-direction. This thus 



generates two droplets and the distance between their the 
centers of mass is denoted by R\2- We have found that 
when i?i2 is bigger than 14.2cr, the two droplets stay as 
separate and stable entities - at least during the time of 
200r and and T = 0.71e/fcs. Otherwise they merge dif- 
fusively as illustrated in Figure [| There appears to be 
no noticeable hydrodynamic velocity field related to the 
coalescence. 
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FIG. 6. The density contours of the droplets during coales- 
cence. The data were time average over periods of 25r. The 
top panels refer to a central cut in the xz plane. The bottom 
panels refer to a plane which is half way between the initially 
separated droplets and later becomes an equatorial plane in 
the fused droplet. 



VI. GRAVITY DRIVEN FALL OF DROPLETS 
AND BUBBLES 

We now return to the problem of the droplet or the 
bubble motion in the pore and consider the case in which 
g is applied only to the atoms of the moving object. We 
ask whether this motion is similar to that of the rigid 
amorphous solid. 

A hydrodynamic prediction |38| for the drag force act- 
ing on a spherical liquid droplet or bubble of radius R 



P 9 p ^ + 3rf 

F d = 2j:Rr] ; — — V 

Tj + r) 



(8) 



where 77 is the viscosity of the surrounding fluid and rj' is 
viscosity of the fluid of which the moving object is made. 
This coincides with the Stokes law (equation [l]) in the 
limit of rf — > oo. In the other limiting case of the gas 
bubble, r/' — > 0, we obtain: 



Determination of the terminal speed results from balanc- 
ing the gravitational pull by the drag force combined with 
the buoyant force and then p3] 



2R 2 g(p'~p)( V + rf) 
377(277 + 377') 



(10) 



where p is the fluid density and p' is the moving body 
density. The sign of Vt agrees with g for droplets but is 
opposite to g for bubbles. In our MD simulations, the 
effects of buoyancy are not incorporated because of the 
imposition of the periodic boundary conditions in the ax- 
ial direction. The corresponding hydrodynamic formula 
for vt corresponds to setting p in equation 10 to zero 
(and then both a droplet and a bubble move along the 
z direction). Due to diffuse density profiles, a more ac- 
curate comparison of MD to the hydrodynamic results is 
obtained by noting that in the stationary state F d should 
be balanced by g times the number of molecules in the 
moving object. The resulting v t will be denoted as the 
hydrodynamic prediction and the values of viscosity are 
as determined by MD in Section II. 

Results of our MD studies of the droplet motion are 
illustrated in Figure [?]. The left panel shows the time 
dependence of the velocity of the center of mass of the 
droplet. The effects of saturation arise after approx- 
imately 100 and 150t for g = 0.01 and 0.05e/er re- 
spectively. The corresponding v t are 0.265 ± 0.024 and 
1.293 ± 0.021cr/r. The droplet is spherical initially and 
then it has a radius of 4.4cr. The right hand panels of 
Figure for the density contours in the stationary state 
indicate that the droplet continues to be nearly spherical 
troughout the motion if g is sufficiently small. On the 
other hand, for g—0.Q5e/(T, there is a noticeable stretch- 
ing along the direction of motion (the axial and sideway 
radii differ by about 2.7 a). 

The hydrodynamic prediction for v t is 0.246ct/t for 
5=0.01 if the initial value of R. This agrees very well 
with the result of MD despite the perturbing effects of the 
pore walls. This suggests that the continuum approach 
works also for motion of microscopic droplets. Note also 
that the values of Vt agree very well with those for the 
similarly built rigid grain - the center of mass moves 
virtually identically. 

The velocity profiles, shown in Figure 8, are quite sim- 
ilar to those shown in Figure 4 except that there is a 
parabolic profiling in the region of the droplet. 



Fd = AnRr/v 



(9) 




FIG. 7. Time evolution of the velocity of mass center of 
the droplet for the two indicated values of g. The two panels 
on the right show the corresponding density contours in the 
stationary state. The thickest of the solid lines correponds to 
a half of the maximum value of the density. 
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FIG. 8. The left panel shows velocity profiles for a droplet 
accelerated by force g with the values (in units of e/cr) as 
indicated. The right hand panel shows velocity profiles for 
the bubble. 

Figure 9 shows that the saturation stage of a moving 
bubble is reached faster than in the case of the droplet 
and that the shape of the bubble gets affected by the 
motion more heavily than the droplet. The bubble does 
not split and moves as a single entity at least for g less 
than 0.025e/o\ The bubble disintegrates for g of order 
0.05e/o\ The terminal velocities are 0.154 ± 0.07 and 
0.385 ± 0.015-r/cr for g = 0.01 and 0.025e/er respectively. 
The hydrodynamic prediction, however, yields v t that is 
about 10 times smaller. This suggests that microscopic 
bubbles do not move hydrodynamically. The hydrody- 
namic behavior might be restored with a much wider 
pore but we had no capacity to test this hypothesis. 



FIG. 9. Same as in Figure 7 but for the bubble. 

Despite the significant shape distortions, the velocity 
profiles (Figure 8) for the bubble look very much like for 
a falling frozen solid of a comparable density - there is 
no parabolicity in the center. 

VII. GRAVITY DRIVEN FLOW WITH 
DROPLETS OR BUBBLES 

We now consider the situation in which the gravita- 
tional force is applied to all molecules in the pore space. 
This is then two-phase gravity driven flow. 

We find that the terminal velocities are now about 
2.7 bigger than the hydrodynamic prediction would yield 
even though the fluid right next to the pore wall is al- 
most pinned. For the droplet, vt is 0.672 ± 0.018 and 
3.272 ± 0.023cr/r for #=0.01 and 0.05 e/a respectively. 
For the bubble, v t is almost the same as when g was ap- 
plied only to the molecules of the bubble. For instance, 
v t is 0.157 ± 0.011cr/r for g = 0.01e/cr. These results to- 
gether indicate that when velocity profiling outside of the 
moving object become significant the drag force deviates 
from the simple form as given in equation 8. 
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FIG. 10. Same as in Figure 7 but for g applied to all atoms 
within the poore space. 




way around. This existence of this effect may depend 
on the geometry and viscosities involved. In each case, 
however, the velocity profile determined in slices which 
are vertically away from the droplet or the bubble agree 
with the single fluid results of Figure 2. 
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FIG. 12. The velocity profiles when g of 0.01 e/a is applied 
to all atoms withing the pore. The top and bottom curves 
refer to the droplet and bubble cases respectively. The data 
points marked by A and B are the same as in the lower panel 
of Figure 2. 

We conclude by noting that the drag force appears to 
have the simple hydrodynamic form if the moving object 
is dense and small compared to the pore diameter. Oth- 
erwise, the velocity field is significantly disturbed by the 
pore and the drag force becomes more complicated. It 
may also be sensitive to the nature of interactions be- 
tween the pore wall and fluid atoms. 
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FIG. 11. Same as in Figure 9 but for g applied to all atoms 
within the pore space. 

The density profiles shown if Figures 10 and 11 are 
close to those shown in Figures 7 and 9 respectively but 
the axial stretching is a bit stronger. 

Figure 12 shows the velocity profiles, in the moving 
frame of reference, when g is applied to all molecules 
within the pore and compares it to the single fluid re- 
sults as replotted from Figure 2. Both the droplet and 
the bubble continue to have velocity profiles which are 
flat in the center. When compared to the situation in 
which g is applied only to the moving object, as in Figure 
8, the the velocity profile of the droplet is significantly 
enhanced. This is consistent with the enhanced values 
of the terminal velocity. On the other hand, the veloci- 
ties around the center of the bubble are slightly reduced. 
It is interesting to note that, in the central region, the 
droplet moves faster than any of the constituting fluids 
at the same value of g. For the bubble, it is the other 
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